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We have extended a simple model of nonlinear diffusive shock acceleration (Berezhko 



& Ellison 199£; Ellison & Berezhko 1999a) to include the injection and acceleration of 



electrons and the production of photons from bremsstrahlung, synchrotron, inverse- 
+3 • Compton, and pion-decay processes. We argue that the results of this model, which 

is simpler to use than more elaborate ones, offer a significant improvement over test- 
particle, power-law spectra which are often used in astrophysical applications of diffusive 
shock acceleration. With an evolutionary supernova remnant (SNR) model to obtain 
shock parameters as functions of ambient interstellar medium parameters and time, we 
predict broad-band continuum photon emission from supernova remnants in general, 
and SN1006 in particular, showing that our results compare well with the more complete 
time-dependent and spherically symmetric nonlinear model of Berezhko, Ksenofontov, 
Sz Petukhov ( |1999a ). We discuss the implications nonlinear shock acceleration has for 



X-ray line emission, and use our model to describe how ambient conditions determine 
the TeV/radio flux ratio, an important parameter for 7-ray observations of radio SNRs. 

Subject headings: acceleration of particles — cosmic rays — supernova remnants - 
radiation mechanisms: non-thermal — gamma-rays: theory — ISM: individual (SN1006) 



1 Universities Space Research Association 



- 2 - 



INTRODUCTION 



More than twenty years of spacecraft observations in the heliosphere have proven that colli- 
sionless shocks can accelerate particles with high efficiency, i.e., 10-50% of the ram energy can go 
into superthermal particles (e.g., Eichler 1981 ; Gosling et al. 1981 ; Ellison et al. 1990| ). A recent 
example comes from Terasawa et al. ( |1999 ) who report on an interplanetary shock, observed by 
GEOTAIL, where the pressures in thermal and accelerated particles, and the magnetic field are in 
approximate equipartition in the foreshock. They observed that the shock transition was smoothed 
by the energetic particle backpressure, giving unambiguous evidence that nonlinear acceleration 
occurred in an interplanetary shock, weak by astrophysical standards (for this exceptionally strong 
interplanetary shock, the sonic Mach number was M$q ~ 5 — 10). Energetic particles exist through- 
out the universe and shocks are commonly associated with them, confirming that shock acceleration 
is important beyond the heliosphere. In fact, shocks are believed to be the main source of Galactic 
cosmic rays, and energetic particles from shocks may produce spectacular events such as 7-ray 
bursts and X-ray and gamma-ray flaring in blazars. The shocks in these objects are expected to 
be much stronger than those in the heliosphere and can only be more efficient and nonlinear. The 
conjecture that collisionless shocks are efficient accelerators is strengthened by results from plasma 
simulations which show efficient shock acceleration consistent with spacecraft observations (e.g., 



Scholer, Trattner, &; Kucharek 1992 ; Giacalone et al. 1997 ). Other indirect evidence comes from 
the implied efficiency of galactic cosmic ray production, radio emission from supernova remnants 
(SNRs) (see Reynolds & Ellison 1992j ), equipartition arguments in AGNs and 7-ray bursts, etc. 



(see Blandford & Eichler 1987 for an early review). There is also clear evidence that shocks can 
produce strong self-generated turbulence. This has long been seen in heliospheric shocks (e.g., Lee 
1982, 1983; Kennel et al. 1984; Baring et al. |1997| ) and there is evidence that it occurs at SNRs as 



well (i.e., Achterberg, Blandford, & Reynolds 1994) 



Despite this compelling evidence for efficient, nonlinear shock acceleration, many astrophysical 
applications of shock acceleration continue to use the test-particle theory proposed more than 20 
years ago by Axford, Leer, & Skadron ( |l977| ), Krymskii fll977|) , Bell ( |l978; ), and Blandford & 
Ostriker ( 1978] ). We believe it is possible that test-particle models are used in situations where 
nonlinear shocks are clearly expected because the numerical, nonlinear results are complicated 
and unwieldy. None of the current techniques modeling nonlinear shocks have simple analytic 
approximations and this makes it difficult to perform parameter searches or compare nonlinear 
results to observations. 

Here, in an attempt to remedy this situation, we present a simple model of diffusive shock 
acceleration, based on more complete studies, which includes the essential physics of nonlinear 
acceleration when the backpressure from energetic ions modifies the shock structure and induces 
spectral curvature. This model is computationally fast and easy-to-use, but complete enough to 
be a valuable tool for interpreting observations. We also investigate some implications of efficient 
acceleration and the associated nonlinear effects on the modeling and interpretation of particle and 
photon observations of SNRs, and include a detailed study of the broad-band continuum emission 
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from the forward shock in SN1006. 



Our basic model is given in Berezhko & Ellison (1999) where only protons are considered. 
Here, we extend the model to include electrons and alpha particles, and calculate the broad-band, 
continuum photon emission from the ion and electron distributions (Ellison & Berezhko 1999a| ), 



For protons, the injection process whereby some fraction of the shock heated plasma becomes 
accelerated to superthermal energies is, of necessity, parameterized but the model allows for the 
investigation of parameters against observations. Unfortunately, the theory of electron injection and 
acceleration in shocks is on a less secure footing than for protons, so we are forced to use additional 
parameters for electron injection. Helium is included without accounting for the enhancement 
effect known to occur for high mass to charge number ions (e.g., Baring et al. 1999j ), but this 
approximation has little effect on the results presented here. 

We claim our nonlinear model, in spite of its approximations and parameters, is far more 
physically meaningful than the test-particle power laws that are still routinely assumed by many 
workers as the outcome from shock acceleration. A crucial property of efficient, nonlinear shock 
acceleration is the interconnection of the entire particle distribution from thermal to the highest 
energies, and the linkage between protons and electrons. Because energy is conserved, a change 
in the production efficiency of the highest energy cosmic rays must impact the thermal properties 
of the shock heated gas and vice versa. If more energy goes into relativistic particles, less is 
available to heat the gas. In contrast, the power laws assumed by test-particle models have no 
connection with the thermal gas, energy conservation does not constrain the normalization of the 
power law, and the spectral index can be changed with no feedback on the thermal plasma. The 
availability of an easy-to-use nonlinear shock acceleration model will have important implications 
for interpreting the broad-band emission from cosmic sources. In particular, the model presented 
here can efficiently explore parameter space to determine which quantities are the most dominant 
in controlling the spectral character and detectability of a remnant's emission in each waveband. 
The density of the ambient interstellar medium and the environmental magnetic field strength are 
the most interesting and critical parameters in this regard. It is clear that adjusting parameters 
to fit one component, say radio, changes the predicted fluxes at all other frequencies and most 
significantly in the X-ray band. Based on constraints imposed by radio and 7-ray observations, 
the simple model can distinguish, as a function of source model and environmental parameters, the 
relative contributions of synchrotron (from TeV electrons) and non-thermal bremsstrahlung to the 
X-ray continuum. Our model predicts the shape of the full, nonthermal electron spectrum and, for 
sources where non-thermal bremsstrahlung is significant, X-ray line models may need to consider 
these non-thermal distributions. 

Conversely, the model possesses the ability to predict TeV fluxes given detections at other 
energies. One of the most perplexing developments of cosmic ray physics is the lack of clear 
detections of pion-decay photons from SNRs. The original predictions of Drury, Aharonian, &; 
Volk ( |1994 ) are well above current upper limits from EGRET and ground-based TeV telescopes. 



It is almost universally assumed that supernova blast waves accelerate cosmic ray ions at least up 
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to the "knee" near 10 eV, and if so, these high energy ions will interact with the ambient gas 
and produce pions which will decay into GeV and TeV photons. The distinctive pion-decay bump 
was expected to be seen with current sensitivities and the lack of detection means there is still no 
direct evidence that SNRs produce cosmic ray ions and is something of a concern for both theorists 
and builders of 7-ray telescopes. Since the detectability depends on source parameters, observers 
critically need selection criteria which reliably predict which SNRs, observed at other frequencies, 
are likely to be bright at 7-ray energies. This is particularly true for ground-based air Cerenkov 
telescopes which must devote weeks or even months observing a single position in the sky to obtain 
good statistics. As we show in our survey below, the model can provide these selection criteria. 
For most of the cases we have studied, supernovas which explode in low-density, low magnetic field 
regions have the largest TeV to radio flux ratios. 



SIMPLE NONLINEAR SHOCK MODEL 



The model of Berezhko & Ellison ( 1999|) synthesizes the essential features of two complimentary 
approaches to nonlinear shock acceleration, namely the semi-analytical diffusion equation method of 



Berezhko et al. (1996) and others (e.g., Kang & Jones 1991, 1995) and the Monte Carlo technique of 
e.g., Ellison et al. dl996|) . The premise of the model of Berezhko & Ellison is that the main aspects of 
time-dependent, nonlinear shock acceleration can be understood in the framework of a plane-wave, 
steady-state assumption with simple approximations for the accelerated particle spectrum. An 
essential element of nonlinear acceleration is particle escape and Ellison & Berezhko Ql999b ) showed 
(and we confirm here) that the dynamic effects of particle escape are essentially the same in a plane- 
wave, steady-state shock as in an evolving, spherical shock. The two very different scenarios agree 
because the process of particle decoupling (i.e., escape) from shocks at high energies produces the 
same shock modification and increased compression ratio in the evolving, spherical shock solution, 
where particle 'escape' is mimicked by energy 'dilution' in the ever expanding upstream volume, 
as it does in the plane-wave approximation, where particles leave at a maximum cutoff energy or 
a free-escape-boundary (FEB). The essential point is that particle escape from strong shocks is a 
fundamental part of the nonlinear acceleration process and is determined primarily by energy and 
momentum conservation, not time-dependence or a particular geometry. The injection efficiency, 
together with the 'shock size' (however it may be determined), set the shock structure and determine 
the overall acceleration efficiency. 

In other words, in a real, finite sized shock, geometry is important for determining the max- 
imum value of momentum, Pniax? 

where the spectrum cuts off, and particle escape comes about 
naturally. In the plane-wave approximation, escape is parameterized with a FEB, but the model 
can still be used to good effect with this restriction. Similarly, since the majority of particles accel- 
erated by an expanding SNR blast wave typically have acceleration times which are much shorter 
than the time scale of the system (in fact, even the highest energy particles expected in SNRs have 



acceleration times a few times shorter than the remnant age, e.g., Berezhko 1996), a steady-state 
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approximation can be used for modeling the particles accelerated by a strongly modified shock. 
We give a quantitative verification of this claim in Ellison & Berezhko (1999b), and demonstrate 
it below in detail for the SNR SN1006. 



2.1. Protons 

The non-linearity of the acceleration process in high Mach number shocks is manifested through 



the feedback of the ions on the spatial profile of the flow velocity (e.g., Drury 1983; Jones & Ellison 



1991), which in turn determines the particle distribution. The accelerated population presses 
on the upstream plasma and slows it. An upstream precursor forms, in which the flow speed 
(in the absence of instabilities) is monotonically decreasing. The net effect that emerges is one 
where the overall compression ratio, r tot , from far upstream to far downstream of the subshock 
discontinuity, exceeds that obtained in the test-particle scenario, while the subshock compression 
ratio, r su b, which is mainly responsible for heating the gas, is less than the test-particle compression. 
This phenomenon was identified by Eichler ( |1984Q , and Ellison & Eichler (|1984| ), and arises for 



two reasons. The most important is that particle escape from strong shocks drains energy and 
pressure which must be compensated for by ramping up the overall compression ratio to conserve 
the fluxes. The second reason is that relativistic particles, with their softer equation of state, 
contribute significantly to the total pressure making the shocked plasma more compressible. Since 
particle diffusion lengths are generally increasing functions of momentum (e.g., Giacalone et al. 
1993 ; Scholer, Kucharek, & Giacalone [2000| ), high momentum particles sample a broader portion 



of the flow velocity profile, and hence experience larger compression ratios than low momentum 
particles. Consequently, higher momentum particles have a flatter power-law index than those at 
lower momenta, thereby dominating the pressure in a nonlinear fashion and producing a concave 
upward spectral curvature which is the trademark of nonlinear shock acceleration. 

It is precisely this spectral concavity that can be compactly approximated by the simple model 



of Berezhko & Ellison (1999). Berezhko & Ellison assume that the accelerated part of the shocked 
particle phase-space distribution, f(p), above a superthermal injection momentum, p m j, can be 
described as a three-component power law: 

{flinj (P/Pini)~ qsub if Pinj < P < m p C , 

Omc \p/(m p c)]~ gint if m p c<p<p 01 , (1) 
"max (p/Poi) 9min A Poi < P < P max 
Here, t/ su b is determined by the subshock compression ratio i.e., 

<?sub = 3r sub /(r sub - 1) , (2) 

ymin 

is given by, 

3.5-0.5r sub 

Qmin = 3.5 + , (3) 

2 r to t - r suh - 1 
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(see Berezhko 1996 for a full discussion of this equation and Malkov 1997, 1999| for an alternative 



derivation) and 

Vint = (q'tp + <7min)/2 , (4) 

where the prime indicates that we use 

, _ ( 72g + l)M s 2 



(72 g - 1)M| + 2 



<4 (5) 



to calculate q' tp . The far upstream sonic Mach number is Mgo, 72 g = 5/3 is the adiabatic index of 
the shocked gas excluding cosmic rays, and the normalization factors are a mc = a- m ^ [pmj/ ("»p c)] 9sub , 
and a max = a mc (m p c/poi) 9int - The full explanation of these terms (including the normalization aj n j 
in terms of input parameters) is given in Berezhko & Ellison (1999), but for here it is only necessary 



to state that the nonlinear spectrum can be determined with four arbitrary parameters: (1) pi n j, 
(2) the rate of proton injection, r7i n j jP , (3) the maximum momentum, p max , where the spectrum 
cuts off, and (4) poi- The results are relatively insensitive to poi which we take as poi = 0.01p max 
if poi > m p c or p m = m p c otherwise, and p- m j and Jfrnj.p are related, leaving essentially two free 
parameters to describe the proton spectrum (here and elsewhere, m p is the proton rest mass). In 
this steady-state description, the proton injection rate, ?7inj,p (Berezhko et al. use the notation 77), is 
equivalent to the fraction of total shocked proton number density, 712, in protons with momentum 
P > Piny The scheme we adopt for determining the scale length that fixes p max is described in 
Section 3.1. 

The thermal portion of the downstream spectrum is taken to be a Maxwell-Boltzmann dis- 
tribution at the shocked temperature and density, i.e., T2 = P2m p /{p2 k), where the downstream 
pressure, P2, and downstream density, pi = po r tot ; are determined by the solution, i.e., they couple 
to the shock dynamics. No extra parameters are required to specify the thermal portion of the pro- 
ton distribution other than to assume it is Maxwell-Boltzmann. In Figure |l] (Fig. 5 from Berezhko 



& Ellison 1999 ), we compare results from the simple model with those of a Monte Carlo shock 
simulation used as a reference. The Monte Carlo simulation calculates the entire spectrum from 
thermal energies upward and has an internally consistent determination of the injection efficiency 
Using this value of f7mj )P in the simple model (the solid triangle is at Pi n j) gives an excellent match 
between the two techniques. The extremely high compression ratios occur because only adiabatic 
heating of the shock precursor is assumed. When Alfven wave heating is included (as it is in all of 
the results presented below except those shown in Figure 0), considerably smaller r to t's result (see 



Berezhko & Ellison 1999| for a discussion of Alfven wave heating) .^] 



2 Helium is included in the calculation of Eqs. (0) in an approximate way by adding a pressure contribution equal 
to nHe/fiH times the energetic proton pressure. This neglects any enhancement effects for heavy ions, but these could 



be incorporated in the simple model with a parameterization such as used by Berezhko & Ksenofontov (1999). Based 



on the work of Baring et al. (1999), we expect the contribution of helium to the pion-decay emission to be no more 



than 50% of the proton contribution. 



- 7- 



2.2. Electrons 

Before describing how electrons are treated in our simple model, we note that electron injection 
is harder to characterize than proton injection for several fundamental reasons. One stems from 
the basic point that the shock structure is determined almost totally by the ions, since they carry 
most of the momentum, so general considerations of momentum and energy conservation constrain 
their behavior. Internally-consistent models of ion injection can be developed without a detailed 
knowledge of the complex plasma processes, and these models have been shown to match spacecraft 
observations of ion acceleration at heliospheric shocks (e.g., Ellison, Mdbius, & Paschmann 1990| ; 



Baring et al. 1997] ) and hybrid simulations (e.g., Giacalone et al. 1997| ). Electrons, on the other 



hand, carry little momentum and don't influence the shock structure substantially so they act basi- 
cally as test particles. Because of this, their injection and acceleration efficiencies (at nonrelativistic 
energies at least) are sensitive to the details of the complex wave-particle plasma interactions, which 
are not well understood. Furthermore, since protons determine the overall shock dynamics, plasma 
simulations wishing to describe electrons self-consistently must include protons self-consistently as 
well, forcing prohibitively large ranges in time and length scales. Compounding the problem, space- 



craft observations of heliospheric shocks have not until very recently (i.e., Terasawa et al. 1999 ) 
detected the injection and acceleration of thermal electrons. This, in contrast to the many obser- 
vations of thermal proton injection and acceleration, has made it impossible to infer the properties 
of electron injection and acceleration from in situ observations, or even, until recently, confirm that 
diffusive shock acceleration of electrons from the thermal background takes place. 

In order to include electrons in the simple model we note that particles with the same upstream 
diffusion length, k/uq (k is the diffusion coefficient and uo is the far upstream shock speed) |] have 
the same acceleration rate. That is, electrons and protons will obtain the same spectral shape if 
they have the same diffusion coefficient K e {p) = K p (p) in a considered momentum range. Since the 
diffusion coefficient k = X(R)v/3 depends not only on the particle rigidity, R, but also on their 
speed, v, which is different for electrons and protons at p <1 GeV/c, one might conclude that, in 
general, the shapes of the accelerated electron and proton spectra are different. However, it is easy 
to show that the electron and proton spectral shapes should, in fact, be very similar at superthermal 
energies. At relativistic energies, E 3> m p c 2 , the interaction with the magnetic field of electrons 
and protons of the same gyroradius (or rigidity) is indistinguishable if we assume any effects from 
helicity are insignificant. Therefore, the shape of the electron spectrum, f e (p), at p 3> m p c is the 
same as the proton spectrum, f p (p). Furthermore, to the accuracy of the simple proton model, 
where we assume that protons with p > m p c have speed c, and those with p < m p c have speed 
v = p/m p , the spectral shapes will be the same for all p > m p c. To the same accuracy, electrons 
with momenta p < m p c have smaller diffusion lengths than protons at p = m p c and they 'feel' only 
part of the shock transition, i.e., the subshock as assumed in the simple model. Therefore, just as 
we assume for the protons, the spectrum of accelerated electrons at p < m p c can be approximated 



3 Everywhere, the subscript '0' implies far upstream values and the subscript '2' implies far downstream values. 
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by the power law f e oc p 9sub , independent of the diffusion coefficient K e (p), as long as k is an 
increasing function of momentum. 

Thus within the accuracy of the simple model, 

fe(p) = K ep f p (p), (6) 

at all superthermal energies, where the value of the numerical factor, K ep , is determined by 
the relation between the electron and proton injection rates, i.e., using Eq. (|l|) we have K ep = 
(oinj,e/oinj,p) [Pinj,e/pinj,p]' ?sub j where the subscripts e and p denote electrons and protons, respec- 
tively. Unfortunately there are no reliable theoretical predictions for this relation so we parame- 
terize it. We remark that if the damping of high frequency waves near thermal energies is strong 
enough to inhibit electron scattering, there will exist a range of suprathermal electron momenta 
where the electron diffusion coefficient and length are larger than the proton ones, i.e., n e > n p 
(e.g., Levinson $99% , |99|). This case is considered in Baring et al. ( |1999j ), where the inefficient 



scattering of electrons causes their diffusion lengths to lengthen and increases their steady-state 
injection efficiencies. If this is the case, the electron spectra in the sub-GeV range may be quite dif- 
ferent from those shown here and these differences can be particularly relevant to radio synchrotron 
and nonthermal bremsstrahlung emission. 

If we make the specific assumption that the scattering mean free path is 

= Vmip ^g,max (^g/^g,max) j (7) 



where r/ m f p is a parameter independent of particle momentum (Baring et al. 199£ , use the notation 
v)> r g = p/ (qB) is the gyroradius in SI units, r giiriax is the gyroradius at the maximum momentum, 
Pmaxj and a is a constant parameter (a = 1 = n m f p is roughly the Bohm limit), then electrons and 
protons which satisfy, 

f?mfp,e Pe V e = ??mfp,p P% V p , (8) 

obtain the same spectral shape. Since we obtain the proton spectrum from the simple model, 
and if we assume that ?7 m fp is the same for electrons and protons, Eq. (^) allows us to determine 
the electron spectrum by finding the momenta dividing the three-component power law, i.e., the 
intervals where q su b, Qintt and g m i n in Eq. (|T|) apply. We note that the above is just for concreteness 
and that, within the approximations of our simple model, Eq. @ holds and we do not have to 
make any specific assumptions concerning the form of the particle diffusion coefficient other than 
that it is a monotonically increasing function of momentum. 

In the top panel of Figure || we show examples of f e (p) (dashed and dotted curves) and f p (p) 
(solid curve) obtained assuming Eq. (0) with a = 1. As just explained, electrons and protons 
have identical superthermal shapes apart from a slight offset near m p c caused by assuming Eq. (|8|) 
rather than Eq. (|6|). 

To set the amplitude of the superthermal electron spectrum, we find it convenient to use 
the electron to proton density ratio at relativistic energies, (e/p) re \. We chose (e/p) re \, rather than 
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some injection parameter, r/; n j je , analogous to f?inj,p> since (e/p) re i can be constrained by observations 
such as those of Galactic cosmic rays, which suggest values around 1-5% (e.g., Miiller & Tang 1987 ; 
Miiller et al. 1995| ) in the 1-10 GeV range. The final parameter needed to completely specify the 
electron spectrum is the shocked electron to proton temperature ratio, T e 2/T p 2, for which there is 
considerable latitude given uncertainties in the shock speed and other parameters for individual 
SNRs that leaves T p 2 ill-determined observationally. Note that detections of X-ray emission (e.g., 
Zimmermann, Trumper, & Yorke 1996| ) can constrain T e 2 if contributions from nonthermal electron 
components can be ruled out. Given T e 2/T p 2 and (e/p) re i, it is straightforward to compute r/mj.ej 
as follows. The number density of particles per unit energy interval, dN/dE, is given by, 

dN a 2 m \ d-P / n 

dE = Wm dE ' (9) 

where dp/dE = (E + tuq c 2 )/(pc 2 ), for tuq = m e ,m p as the case may be. Here and elsewhere, E 
is the kinetic energy. Exploiting the fact that the electron and proton phase space distributions 
have essentially the same shape we can write the ratio at a given nonrelativistic energy, E e = E p = 
Enr <^m e c 2 , 



Vinj,e _ (dN/dE) c 



r/inj, P (dN/dE) 



amj, C m e p e (Pe/Pm),e) 9s " b _ «inj,e / m e \ / T c ; X '" >> " ' 



l in],p rr >'pPp (Pp/Pinj,p) 9sub °inj,p \ m p J \^p2 



fc2 



(10) 



since Pi n j,e oc (m e T e 2) 1//2 , p e oc (rUeEe) 1 / 2 and likewise for protons. This can be inverted to yield 
ainj,e/oinj,p- At identical relativistic energies, E e = E p = E TC \ ^> m p c 2 , since poi is identical for 
protons and electrons, we have, 



(dN/dE) r 



Pe "max,e iPp/POl) "max,e ^inj,e ( Tie 

PpO maXiP (p e /P0l)~ qm[a a max,p ainj,p \m p 



<7sub/2 

, (11) 



since p p = p e at relativistic energies. The last equality is obtained using the definition of a max given 
just after Eq. (|5|). Eqs. ( |l0| ) and (^) define the relationship between fftnj.e, T e 2/T P 2, and {e/p) Te \, 
with the explicit dependence on T e 2/T p 2 disappearing when ai n j >e is specified. 



Comparison of Eqs. ( |T0[ ) and ( jTl| ) shows that the number density ratio decreases by a factor of 
(m e /m p )( qsuh ~ 3 ^ 2 in going from fully nonrelativistic to fully relativistic energies. For strong shocks, 
r sub ~ 3 and q su ^ ~ 4.5, so producing an electron to proton ratio of (e/p) re \ ~ 0.01 at relativistic 
energies, as observed for galactic CRs, implies that the electron injection rate at nonrelativistic 
energies must be approximately equal to the proton rate, i.e., r/i n j,e — f/inj.p- This conclusion holds 
whether the thermal electrons and protons are in equipartition or not. However, it is highly likely 
that the injection energy scales strongly with the shocked temperature, implying that T e i ~ T p 2, 
i.e., it is unlikely that (e/p) ve \ ~ 0.01 can be obtained unless the shocked electron and proton 
temperatures are nearly equal a result suggested by some thermal X-ray signatures of SNRs (e.g., 
Zimmermann, Trumper, & Yorke |199C| ). 



4 The work presented here implicitly assumes that all of the accelerated protons and electrons originate from the 
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For the example shown in Figure [| we have chosen T e 2/T p 2 = 1, and used (e/p) m \ = 0.05 
(dashed line) and (e/p) ie \ = 0.01 (dotted line). The spectra also include an exponential cutoff (not 
included in the initial Berezhko & Ellison [1999| work) so that f(p) given by Eqs. (|l]) becomes (see 
Berezhko & Krymsky 1996 and Berezhko et al. 1996[ ): 



f(p) -> f(p) exp 



a \p 



P 



(12) 



where the a is the same as in Eqn. (0). 



The bottom panel of Figure || shows the differential flux, dJ/dE = [v / (4ir)]dN/ dE, in particles 
per cm 2 per sec per ster per MeV. For {e/p) ve \ = 0.05, ??mj,e = 0.012, while for (e/p) re \ = 0.01, 
r^inj e = 2.5 x 10~ 3 , if we assume that the electron and proton injection energies are the same. 

Finally, we define the acceleration efficiency, e re i, as the fraction of total incoming energy flux 
that goes into particles with momentum p > m p c, i.e., 



Crel 



[72/ (72 - 1)] -Prel,2(> ™p C) U 2 + F E 

(l/2)p «g + (5/2) P u 



(13) 



where -P re i,2(> m p c ) is the downstream pressure in particles of momentum m p c or greater, u is the 
bulk flow speed in the shock frame, 72 is the effective downstream ratio of specific heats including 
the effects of relativistic particle pressure, and we have taken the ratio of specific heats for the 
unshocked plasma to be 5/3. The escaping energy flux, Fe, is included in this definition since 
these are accelerated particles that only leave the system after obtaining an energy E ~ E max . The 
efficiency for selected examples are listed in Tables 1 and 2 under 'Output values'. 



2.3. Photons 

The photon emission processes considered in this paper are those discussed at length in Baring 
et al. ( 1999| ): synchrotron emission that pertains to the radio to X-ray wavebands, inverse Compton 



scattering of cosmic microwave background radiation that contributes to the gamma-ray signal, X- 
r ay /gamma-ray bremsstrahlung emission of thermal and energetic electrons interacting with both 
ambient electrons and protons,^ and pion decay gamma-rays spawned by p-p collisions involving 
the accelerated cosmic rays. As in Baring et al. (|1999| ), we explicitly include the photon production 



thermal background. In fact, superthermal seed populations exist, such as the Galactic cosmic rays and, in young 
SNRs, positrons and electrons from radioactive decay of freshly synthesized material will be accelerated if they 
encounter local shocks. For simplicity, however, we assume that the contributions from these sources to the freshly 
accelerated thermal material is small, noting that the ambient ISM density exceeds that of cosmic rays by several 
orders of magnitude. 

5 We note that, to correct some typographical errors and improve the accuracy of an approximation to the e-e 



bremsstrahlung cross-section, Eqs. (A2) and (A4) in Baring et al. (1999) should be modified such that all 7 e 's in 
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from He-p and He-He collisions. While Sturner et al. (1997) were the first to examine the radiative 
signatures of all of these components, Baring et al. (1999) were the first to treat them all in 
the context of a detailed model of nonlinear shock acceleration. Generally, infrared radiation 
fields peculiar to individual remnants contribute only a minority of the inverse Compton signal 
(e.g., Gaisser, Protheroe, & Stanev 1998 ; but see de Jager & Mastichiadis' |1997 discussion of the 



environment of W44). To date, and also in this paper, synchrotron self-Compton contributions are 
neglected, being important only in strongly magnetized (i.e., B ;> 100/xG) remnants such as Cas A. 
Furthermore, inverse bremsstrahlung from fast ions contributes insignificantly in general (Baring, 
Jones, & Ellison |2000 ), and is accordingly neglected here. 



Losses from synchrotron and inverse-Compton emission are included and combined such that 
electron spectra cutoff at an energy 
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(15) 



where is the shock speed and -B CDr = 3.23 pG is the equivalent magnetic field of the cosmic 
background radiation and accounts for inverse-Compton losses. This expression differs slightly from 
that in Baring et al. because we explicitly allow for the downstream field B 2 to be different from 
Bq. Also, in order to be consistent with the work of Berezhko, Ksenofontov, & Petukhov( 1999a| ) , 
we take the perpendicular component of the downstream magnetic field to the line of sight to be 
B%± = O.5B2 in calculating the synchrotron emission. 



3. MODELING SNR PHOTON EMISSION 

Instead of arbitrarily choosing the shock speed and maximum acceleration energy, we use a 
model of SNR evolution by Truelove & McKee ( 1999| ) that continuously maps between the free 



expansion and Sedov phases to give the shock parameters as a function of explosion energy, E sn , 
ejecta mass, M e j, and remnant age, t S nr- For this work, we only consider uniform distributions for 
the ambient ISM mass, appropriate to Type la progenitors. We do allow for a power-law density 
profile for the ejecta, as in Chevalier ( |1982j ), i.e., 



p oc i- n 3 r (r/i snr )- n , (16) 



Eq. (A2) are replaced by 7 e — 1, and Eq. (A4) there is written as 



10 ( 7 c - 1) 1/5 
3 7c + 1 

for e 7 where A(ej, j e ) is positive; otherwise A{e 1 , 7 C ) is set to zero 



^^^-y^Ti-i;) (14) 
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where p is density, r is radius, and n is a constant, but only consider the forward shock. Examples 
with forward and reverse shocks are given in Ellison ( pOOOj ). We caution that care must be used 
in applying our acceleration model and associated predictions of broad-band emission to Type II 
supernovas, which may explode in environments that are far from uniform, being strongly modified 
by the stellar wind generated by their massive progenitors. 

In Figure |3|a we show the instantaneous particle spectra produced by the forward shocks in 
two examples which differ widely in ambient parameters. We have a high density, high ejecta 
mass, high S-field example ('FT with heavy-weight curves) with proton number density n P Q = 10 
cm -3 , M e j = 5Mq, and Bq = 30 x 10~ 6 G (more typical of Type II SNas), and a low density, low 
ejecta mass, low S-field example ('L' with light-weight curves) with n p o = 0.01 cm -3 , M e j = IMq, 
and Bq = 3 x 10 -6 G (more typical of Type la SNae). A full listing of the parameters is given 
in Table 1 under columns marked 'High' and 'Low' but all other input parameters are the same 
between these two models, namely: E sn = 10 51 erg, t snr = 1000 yr, T e 2/T p 2 = 1, (e/p) re \ = 0.03, 
and ?7inj,p = 1 x 10~ 3 . Figure ^b shows photon fluxes at the Earth assuming a distance, D snr = 1 
kpc and an emission volume, V em i s = 1 pc 3 . All output parameters are given in Table 1. In all of 
the examples in this paper we assume the plasma is fully ionized and the ambient helium number 
density to be 0.1 n p Q. We further assume that the unshocked helium temperature per nucleon equals 
the proton temperature, T p q (this is typically observed in the heliosphere), and that the unshocked 
electron temperature equals T p q. 

There are some striking differences between the photon spectra for these parameter sets which 
represent a reasonable spread for various types of SNRs. First, the three order of magnitude 
difference in n P Q and one order difference in Bq produce much stronger emission, per unit emission 
volume, for the high values. Next, for high values, the TeV emission is dominated by pion-decay 
with inverse-Compton being completely insignificant (and not shown in Figure ||b), while for low 
values, the TeV emission is dominated by inverse-Compton. The inverse-Compton from the low 
value example also yields the highest photon energies. At hard X-ray energies (10-100 keV), the 
emission from the high density example has a strong nonthermal bremsstrahlung component, while 
for the low density example, the X-rays < 100 keV are almost totally from synchrotron emission 
from high energy electrons. Deducing the shifting dominance of the various emission processes as 
parameters vary is an important goal of this work. 

3.1. Parameter Survey 

As indicated by the previous examples, nonlinear shock acceleration is complicated and the 
solutions we describe have a large number of parameters, even in the homogeneous environment 
we assume. The input parameters are: E sa , M e j, t snr , n P Q, Bq, n, D snr , and V^ m i s . Added to these 
are the parameters intrinsic to the model: r/mfp, 7/inj,p> T e 2/T p 2, and (e/p) re i, plus what we call 
'general model' assumptions. These include a, pq±, equation (|i~2|), and the assumptions of ambient 
and ejecta mass distributions. The maximum momenta where the spectra cut off are also model 
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parameters, but they are determined from the SNR model (fjmfp becomes important here) as was 
described in Baring et al. ( |1999| ). That is, in determining the maximum momenta of protons, we 
first assume equation (|?]) and then set the diffusive shock acceleration time, i acc , equal to i snr , or 
set the upstream diffusion length equal to 1/4 of the shock radius, whichever produces the lowest 
Pmax-0 For electrons, p max is equal to that of the protons or to the value determined from combined 



synchrotron and inverse-Compton losses (see Baring et al. 199S , for details), whichever is less 



3.1.1. Emission volume 

For a uniform environment, we can estimate the emission volume in the following way. The 
contact discontinuity separates the swept up ISM material from the ejecta and acts as the piston 
for the forward shock. During the free expansion and early Sedov phases, the swept up material 
forms a nearly uniform, dense shell behind the forward shock such that, 

R^kPO ~ Pshell (#sk ~ -^piston) i ( 17 ) 

where i? p i s ton is the radius of the contact discontinuity, po is the density of the background ISM, 
and Psheii is the density in the shell. In the approximation that the shell has uniform density, 
Psheii = ?*totPo> i- e -> it nas the shocked density, p2- For large compression ratios, the thickness of 
the shell, L = — -R p i s ton <C R s k, and we can write, 

J 1 / -Rpiston Y = i f R s k~ L \ 3 ^ 3L ^ 



not V R s k J \ R s k J Rsk 

During later phases, the contact discontinuity will be far from the forward shock but we can then 
interpret L as the thickness of the dense shell behind the forward shock where most of the emission 
occurs independent of the position of the contact discontinuity. With these approximations, the 
emission volume can be estimated as 

Kmi 8 ~(47r/3)i4/r t ot , (19) 
so that the emission volume generally is considerably less than the total remnant volume. 



3.1.2. Trends and Key Parameters 

In Figures^, [|, and [6| we plot the total photon emission varying a single parameter as indicated. 
All of these results are normalized to D snr = 1 kpc and V^ m j s as given in Eqn. (|i~9|). Varying r/i n j,p, 



6 In determining p max at time t snl , we make the approximation that the shock speed is constant at the value it 
has at t sm . In fact, an accurate determination of p max in an evolving remnant requires a more complete model than 



we use here, i.e., Berezhko (1996), which keeps track of the history of particles, adiabatic losses, and the numbers 
of particles accelerated at a given epoch. Despite this, our values will not differ by large factors from more realistic 
ones. 



-14- 



T e 2/T p 2, or (e/p) re \ (Fig||) produces relatively modest changes in the photon emission, particularly 
in the radio to X-ray and gamma-ray bands. An exception occurs for 7?i n j,p = 10~ 5 (solid curve, 
Fig. ^|a). This injection rate is low enough that a test-particle result occurs even though the Mach 
number is quite high (Mso — 130), i.e., the energy placed in superthermal particles is not enough to 
smooth the shock. A test-particle solution is different in two main ways from the nonlinear cases: 
the overall compression ratio, r tot , is approximately four which is considerably less than for the 
higher injection rate cases, producing the steeper spectrum at the highest energies, and the shocked 
temperature is considerably higher than the nonlinear temperatures (as indicated by the position of 
the UV/X-ray thermal bremsstrahlung peak in panel f|a). While the gamma-ray flux drops in the 
test-particle regime, the concurrence of curves for higher ??inj,p underlines the regulatory effect (e.g., 
Eichler 1984; Ellison & Eichler 1984; see also the review of Jones and Ellison |1991| ) of nonlinear 



acceleration. 

Fugure |]b reveals that by itself, the value of the parameter T e 2/T p 2 is of little significance. This 
is because of its coupling to Tftnj.p, and (e/p) re \, as described in Section |2.2| , so that by fixing ??inj,p) and 
(e/p) re \, the non-thermal X-ray/soft gamma-ray bremsstrahlung flux and the hard gamma-ray pion 
decay and inverse Compton contributions remain approximately fixed. The proton temperature is 
determined by the shock dynamics solution. 

Changes in the electron to proton ratio at relativistic energies, (e/p) re \ (Fig. also produce 
important effects in the superthermal bremsstrahlung range, since this component depends on 
the normalization of the non-thermal electron distribution. As {e/p) Te \ increases, the super-GeV 
gamma-ray band shifts from being pion-decay dominated to having nearly equal bremsstrahlung, 
inverse-Compton, and pion-decay components. We show below examples with a more tenuous ISM 
where inverse-Compton dominates the gamma-rays so that the flux in this energy range becomes 
quite sensitive to the choice of (e/p) re \. 

The ambient density is expected to have a large natural range (Fig. ||a) and strongly influences 
the emitted flux. Low values of n p o give the highest maximum energy with the >TeV flux dominated 
by inverse-Compton. As n p Q increases, the pion-decay component becomes more pronounced in the 
100 MeV to TeV range and bremsstrahlung becomes dominant in the X-ray band. X-ray emission is 
dominated by synchrotron at low n p o, and the radio emission is relatively insensitive to n p o with a 
slight flattening occurring at higher densities. Naively, one expects that the synchrotron and inverse 
Compton components should scale linearly with density increases, while the bremsstrahlung and 
pion decay contributions should be proportional to n p0 . However, the global spectral properties 
are an accumulation of effects caused by the complexity of the nonlinear acceleration mechanism 
and the evolution of the SNR as given by the Truelove & McKee parameterization. As n p Q rises, 
the expanding supernova sweeps up its ejecta mass sooner, and therefore decelerates on shorter 
timescales, thereby reducing both the radius and volume of a remnant of given age, and lowering 
the shock speed, the compression ratio, and the downstream proton temperature T p 2 [e.g., see 



Eq. (7) of Baring et al. 199£ ; and Table 2] . The bottom line is that the normalization of the overall 



distribution and the relative contributions of the various components, are not simple functions of 
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density and can vary strongly in narrow energy bands. The full set of parameters for Figure |H| is 
given in Table 2. 

Some remarkable features in Figure ||a (where B = 3/xG for all four examples) include (i) for 
n p o > 1 cm~ 3 (i.e., the dotted curve, see also Figure ||), the radio to optical/UV emission is of 
thermal bremsstrahlung rather than synchrotron origin, (ii) the highest energy emission occurs for 
the lowest n P Q and is dominated by inverse-Compton, and (iii) the spectral shape in the X-ray band 
(1-100 keV) is highly sensitive to the density (thereby providing powerful observational diagnostics) 
and can have a strong nonthermal bremsstrahlung component at high n p Q. 

Variations caused by adjusting the ambient magnetic field are illustrated in Fig. ||b (where 
n p o = 1 cm -3 for all four examples). The principal property immediately apparent is the anti- 
correlation between radio and GeV-TeV fluxes. The higher the value of Bq, the brighter the 
radio synchrotron [scaling as ~ B 3 / 2 : e.g., see Eq. (6.36) of Rybicki & Lightman |1979(| , which is 



the effect usually incorporated in test-particle SNR models (e.g. Mastichiadis & de Jager 1996 ). 
In addition, however, the hard gamma-ray emission becomes fainter with increasing Bq. This 
property arises mostly because of the influence of the magnetic field on the shock dynamics and 
total compression ratio. The Alfvenic Mach number drops and Alfven wave heating in the precursor 
becomes stronger as Bq increases causing r to t to decrease. This weakening of the shock steepens 
the particle distributions and the overall photon spectrum, with the net effect of a reduction of the 
gamma-ray flux together with a change in its spectral index.0 The emission from radio through 
the X-ray band is dominated by synchrotron emission from energetic electrons for Bq > 10 ^G, 
with bremsstrahlung dominating X-rays for Bq < 3 fiG. Another consequence of a rise in Bq is that 
the maximum ion energy increases while the maximum electron energy declines due to synchrotron 
losses. This causes pion-decay to dominate the gamma-ray band and the UV/soft X-ray synchrotron 



peak becomes approximately independent of Bq, a well-known property (e.g., de Jager et al. 1996 



Reynolds 1996) ) of cooling-limited synchrotron radiation. Note again that we omit the regime 



Bq > lOO^G where the Alfvenic Mach number is as low as a few and where second-order Fermi 
acceleration is important (see Bykov et al. 2000| and references therein), an effect beyond the scope 
of the present work. 

The plots in Fig. [| include the approximate integral flux sensitivities (adjusted for the E 2 times 
differential flux representation) of the proposed Veritas experiment (Weekes et al. 1999| ) and the 



planned GLAST mission (Gehrels & Michelson 1999| ). These represent the probability that each 



of these experiments will detect photons above a given energy, and clearly indicate that GLAST 
and Veritas will make significant observational progress in the 100 MeV-100 GeV and 100 GeV-10 
TeV bands, if SNRs emit according to the predictions here. CELESTE provides somewhat better 
sensitivity in the limited 50-70 GeV range, but we chose the Veritas sensitivity as a benchmark for 



7 The spectrum is further steepened because the Alfven waves which do the scattering are assumed to move through 
the upstream plasma such that they move with the Alfven speed away from the shock when viewed from the local 
plasma frame. 
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atmospheric Cerenkov telescopes (ACTs). Several new TeV programs are in the works, including 
MAGIC (Lorenz |l~997j) in the Canary islands, and HESS (e.g. Kohnle et al. |l999|) and CANGAROO- 



III in the southern hemisphere; their sensitivities differ only modestly in energy range and flux from 
the ones depicted for Veritas. Because of the "saturation" of gamma-ray spectra at low densities, 
due to the insensitivity of inverse Compton emission to ISM density, for reasonable ISM fields these 
experiments should have no difficulty detecting remnants throughout the Sedov epoch and prior 
to the radiative phase. The only way this cannot materialize is for shocks at remnant shells not 
to be efficiently accelerating particles to energies beyond around 10 GeV, i.e., contrary to general 
expectations. Therefore, we anticipate positive detections in a number of SNRs in and out of 
the galactic plane, and even the most pessimistic case of no detections will lead to an improved 
understanding of particle acceleration. 

In particular, remnants in moderately dense environs (i.e. ^ 1 cm -3 ) will be ideal candidates 
for GLAST to search for the 7r° bump emission so clearly evident in the higher density cases in 
Fig. |5]a. A clearer indication of this can be found in Baring ( pOOOj ) , where these model spectra are 
compared with the simulated differential flux sensitivity expected for GLAST, a truer indicator of 
the ability of the instrument to perform spectral measurements; Baring ( 2000| ) observes that the 7r° 



bump will be a clearly discernible spectral structure when ISM densities exceed around 0.3 cm -3 . 
Confirmation of the existence of accelerated ions that spawn such a gamma-ray signature is an 
unambiguous signature of cosmic ray production. We note that the EGRET sensitivity over the 
100 MeV - 10 GeV band is in the 10 _5 -10 MeV cm -2 sec -1 range, so that only remnants with 
ambient densities exceeding around 0.3 cm~ 3 would be possibly detectable by EGRET. Moreover, 
EGRET would only be sensitive to pion decay bumps when the environmental densities are quite 
large, well above 10 cm -3 . 

In Fig. ^ we show variations in the SN explosion energy and the age of the remnant. Ad- 
justing E sn clearly produces little change in the broad-band spectral shape, merely modifying the 
normalization due to changes in the radius and volume of the remnant. On the other hand, at early 
times (i.e., i snr ^ 300yr), the SNR age has a strong effect on the flux, mainly because the remnant 
volume is small prior to the Sedov phase. During these early times, which consists mostly of the 
free expansion phase, the shock is extremely nonlinear in nature, obtaining a very large compres- 
sion ratio (r to t — 25 when t SDI = 30 yr and r to t — 24 when i snr = 100 yr),0 and the gamma-ray 
band is dominated by pion-decay. The shock speed is high during this time and the downstream 
temperature, while very much less than the R-H value (T2 ~ 3 x 10 K for the nonlinear shock at 
i snr = 100 yr versus T2 — 1.4 x 10 9 K for the test-particle shock!), is also high. The transition into 
the Sedov phase at t snr ~ 300 yr brings about a steady decline in T p 2 along with a gamma-ray flux 
which is almost independent of SNR age, a property present in the more complex models of Dorfi 



3 Note that high Mach number unmodified shocks (i.e., those with rtot — 4) are possible as described in general 



by Berezhko & Ellison (1999) and specifically for SN1006 by Berezhko, Ksenofontov, & Petukhov (1999a) if injection 



rates lower than that assumed here are used. 
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(|l99lD , Drury, Aharonian, & Volk ( |l994|) , Baring et al.( |l999| ), and Berezhko & Volk (|l997| ). The 



origin of this gamma-ray insensitivity to i snr prior to a remnant's radiative phase is an approximate 
compensation (Baring 2000) between the volume that scales as isnr (radius oc isnr) in the Sedov 



—3/5 

phase, and the shock speed (and therefore also the square root of T p 2) that declines as i S nr so 
that a roughly E~ 2 particle distribution function has a normalization that scales as T p 2 oc i S nr • 
The contrast between the free expansion and Sedov phases is underlined by the change in shock 
parameters: at t S nr = 30 yr (solid curve) V sk = 1.4xl0 4 km s _1 , M so = 850, M A0 = 2.4xl0 3 , 
and rtot = 25, while at i S nr = 10 yr (dotted curve), V; k = 490 km s _1 , M so = 31, M A0 = 90, and 
r tot = 8.5. 

In Figure we show, as a function of Bq and n p o, the following ratio: 

Rphot = F V {E = 10-UMeV) ' (20) 

that is, the integral photon flux above Eq normalized to .Eo over the radio flux (F v has units of 
cm~ 2 -s~ and we set Eq = 500 GeV in Fig. [?]). This ratio is critical for gamma-ray searches of radio 
SNRs and the most important result is that, for n p Q < 1 cm~ 3 , SNRs in high Bq environments 
have the lowest .Rphot- For n p o held fixed at 0.01 cm~ 3 , i? pno t decreases by more than two orders 
of magnitude as Bq increases from 1 to 100 fiG. However, the absolute TeV flux depends strongly 
on the density as shown in Figure |5|a. 

Another critical factor for TeV observations involves the relative importance of inverse-Compton 
versus pion-decay emission and the impact this has on our ideas for the origin of cosmic rays. In 
Figure || we show the inverse-Compton and pion-decay components for the parameters marked 
with dots in Figure [7]. The variation in total absolute flux for the three cases is relatively small 
but there is a shift from inverse-Compton being dominant at TeV energies when n P Q = 0.01 cm -3 
to pion-decay being dominant when n p Q = 10 cm -3 . In the high Bq case, pion-decay emission 
extends to higher energies than the inverse-Compton because the electron spectrum is cutoff from 
synchrotron and inverse-Compton losses. Hence, atmospheric Cerenkov experiments may provide 
detections of pion decay emission in sources with high magnetic fields. This may be the case in 
the recent marginal detection (Volk et al. 2000) ) of Cas A by HEGRA. Yet the spectral shape in 



the TeV band bears no signature peculiar to pion-decay emission, but rather is a marker of the un- 
derlying particle distributions. Hence, the most powerful diagnostic the sensitive TeV experiments 
will provide is the determination of the maximum energy of emission, thereby constraining the ISM 
density, magnetic field and the electron to proton ratio. 



3.2. SNR SN1006 



For a concrete example, we apply our model to SN1006, the first remnant to have a detection 
of shell-related TeV gamma-ray emission reported, by the CANGAROO collaboration (Tanimori 
et al. 1998). We constrain our parameters to match those used by Berezhko, Ksenofontov, & 
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Petukhov (1999a) as closely as possible and discuss the comparison with that model below. A 



full listing of the parameters is given in Table 1 under the column labeled SN1006nl> where the 



'NL' subscript contrasts this nonlinear result with the test-particle one discussed in Section 3.3 
below. The parameters used by Berezhko, Ksenofontov, & Petukhov are listed in the column 
labeled BKP99. In particular, we use an unshocked proton density n p o = 0.1 cm" 3 , SN energy 
E sn = 10 51 erg, ejecta mass M e j = 1.4Mq (i.e., values similar to those given by Truelove & McKee 



1999 and/or Laming et al. 1996 ), and an ejecta density profile characterized by n = 7 (i.e., Eq. 16). 
For a remnant age t snr = 990 yr, these values yield ~ 3750 km s _1 and ~ 7 pc in the 
Truelove and McKee model. If we take the observed proper motion of 0".30 ± 0".04 yr" 1 (Long et 



al. 1988| ), this speed yields a distance of 2.3 kpc, and using the observed angular size of 30', this 



implies an angular radius of 10 pc, within 50% of the 7 pc given above. 

Figure |^ shows the continuum component photon spectra obtained from these parameters, 
along with the sum of the components (heavy solid line). In obtaining the normalization for this 



fit, the emission volume was taken to be ~ 195 pc , i.e., the value given by Eq. (19) and the distance 



to the source was assumed to be D snr = 1.7 kpc, commensurate with that deduced from optical 



observations (see Green 1998, and references therein). If we had taken L> snr = 2.3 kpc, the required 
emission volume would be correspondingly larger to match the observed flux. It's clear that a 
reasonable fit to the observed spectra is produced considering that we used the same parameters 
as Berezhko, Ksenofontov, & Petukhov .[] Note that we have not included free-free absorption in 
our models. 

Our results indicate that the TeV emission in SN1006 is dominated by the inverse-Compton 
component, in general agreement with the conclusions of previous work that invoke test-particle 
power-laws, i.e., Mastichiadis & de Jager (1996) and Tanimori et al. ( |199§| ). However, there is 



a pion-decay component at TeV energies and small changes in the parameters can make it more 
pronounced than shown in Fig. B as indicated in the results of Berezhko, Ksenofontov, & Petukhov 
(1999a, b) discussed below. Also in agreement with previous work (e.g., Koyama et al. 1995 



Reynolds [i996| ), theX -ray emission is dominated by synchrotron emission from super- TeV electrons. 



However, bremsstrahlung does contribute some at hard X-ray energies and again, this contribution 
is sensitive to the input parameters. The result shown in Figure [9| assumes T e 2/T P 2 = 1. AXAF and 
XMM will probe the X-ray emission at and below the ASCA data range depicted in Figure |, and 
the INTEGRAL experiment will be sensitive to the hard X-ray band below 200 keV. Together, these 
spacecraft should be able to clearly differentiate the synchrotron and bremsstrahlung components 
in the X-ray range and, therefore, provide constraints for the entire broad-band emission. 



9 We have used a = 0.6 for this fit. Berezhko, Ksenofontov, & Petukhov ( f 999a| ) and Ammosov et al. (1994) found 
that a similar value (a = 0.5) was required for a good fit and Reynolds (1996) gave a detailed explanation of why 
such a broadening might occur for the electron spectrum. 
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3.2.1. Comparison with the Kinetic NL Model of Berezhko et al. 



As mentioned above, Berezhko, Ksenofontov, & Petukhov ( |1999a , 1999b| ) have presented results 



for SN1006 using a kinetic, nonlinear model of diffusive cosmic ray acceleration by a spherically 
symmetric, expanding SNR. This model self-consistently solves the cosmic ray transport equations 



together with the gas dynamic equations and is based on the work of Berezhko et al. ( 1996 ) and 



Berezhko Sz Volk (1997). It is more complete than the steady-state, plane-shock model presented 
here and keeps track of the accelerated particles as the SNR evolves and yields the integrated emis- 
sion over the entire remnant as a function of time including adiabatic losses. It also considers the 
diffusive properties of the particles as a function of energy so that low energy electrons responsible 
for the radio emission occupy a thin region behind the forward shock, while the highest energy 
particles fill the entire volume behind the shock. Despite these substantial differences, both models 
produce good fits to the broad-band continuum with virtually identical input parameters. The 
Berezhko, Ksenofontov, &; Petukhov parameters are listed in Table 1 under the heading 'BKP99.' 

The only difference in the input parameters is the small difference in a, all other input param- 
eters are the same. Note that the proton injection rate we use (?7mj,p = 2.1 xlO -4 ) is equivalent to 
the value (?7inj,p = 5xl0 -4 ) used by Berezhko, Ksenofontov, & Petukhov, since injection is treated 
somewhat differently in the two models and these two values give the same injection efficiency. The 



output parameters, as indicated in Table 1, differ only slightly and in Figure 10 we compare the 
shapes of the proton distributions for the two models. The most obvious difference is the lack of 
a thermal peak in the Berezhko et al. result (solid curve). They inject particles at a superthermal 
momentum (i.e., at p- m j as given in Eqn. |l]) and do not model thermal particles. However, they 
could have included a shock heated thermal distribution as we do here, in which case it would 
have been very similar to the dashed curve since both models have similar pre-shock and injection 
conditions. The good correspondence in shape over the entire superthermal range occurs despite 
the fact that the plane-shock, simple model is steady-state and uses the shock conditions at the 
SNR age (as determined by the Truelove & McKee parameterization) to determine the emission 
at that instant. In contrast, the Berezhko et al. model is fully time-dependent and integrates the 
emission over the remnant. Notably, the spectral curvature typical of nonlinear shock acceleration 
survives the volume integration in the Berezhko et al. result. 

The two models do differ somewhat in the fraction of TeV flux contributed by pion-decay. 
Berezhko, Ksenofontov, & Petukhov find that pion-decay and inverse-Compton contribute about 
equally while here, inverse-Compton is clearly dominant. While these differences and other effects 
will be considered with more detailed comparisons in future work, we believe this example gives 
credence to our claim that the simple model is getting the basic physics correct. Furthermore, 
the added ability to provide rapid, broad-band fitting is very important because it is simply not 
possible, because of computing limitations, to perform a comprehensive parameter search with 
complicated codes such as the Monte Carlo technique used in Baring et al. ( 1999j ) in a timely 



manner. Our model does permit such a search, and we anticipate it can be used to significant 
advantage in astronomical data analysis packages, providing the capability, not previously available, 
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for incorporating the principal nonlinear effects of shock acceleration. 



3.3. Test-particle vs. nonlinear results 

As we mentioned above, despite considerable evidence to the contrary, many astrophysical 
applications of diffusive shock acceleration still assume that the acceleration is inefficient and pro- 
duces test-particle power laws.[^] A high efficiency for producing cosmic rays will have far-reaching 
consequences for a number of applications because values for the shocked temperature and density 
will differ from the test-particle Rankine-Hugoniot values. This will be particularly important for 
X-ray line emission models. To illustrate how the predictions of diffusive shock acceleration vary 



depending on whether the test-particle or nonlinear assumption is made, we show in Figure |11 
(light-weight curves) the phase-space electron and proton distributions that produced the emission 
shown in Figure [| These are compared to test-particle spectra (heavy-weight curves) calculated 
for parameters which are the same except the proton injection efficiency is set to r)m),p = 1 x 1CP 5 
in the test-particle case. 

For the test-particle case (column labeled SN1006tp in Table 1): r to t = 4.1, r su b = 4.0, 
and T2 — 1.7 x 10 8 K, while for the nonlinear case: r tot = 7.2, r su b = 3.6, and T2 — 4.7 xlO 7 
K. In the nonlinear case, the shocked density is ~ 75% greater and the shocked temperature is 
more than a factor of 3 less than the test-particle result. Furthermore, the nonlinear electron 



spectrum in Figure 11 shows a stronger nonthermal tail on the thermal peak (comparison between 
the dotted and dot-dashed curves in Figure Q gives a much more extreme example). Accurate 
electron temperatures and knowledge of the presence of significant nonthermal electron distributions 
are of prime importance for modeling X-ray line emission from shock heated gas (e.g., Hamilton 
et al. 1983; ). Line models depend strongly on the assumed electron temperature and density and 



may change substantially if cosmic rays are accelerated with efficiencies typical of the example just 



given. Dorfi & Bohringer (1993) (see also Dorfi 1994) calculate the evolution of SNRs including 



cosmic ray production and find examples where ;> 50% of the explosion energy is transferred to 
cosmic rays, reducing the X-ray emission by a factor of 20 from test-particle predictions. Efficiencies 
nearly this high are found by Berezhko, Ksenofontov, & Petukhov ( |1999a , 1999b| ) and are typical 
of those shown here in Table 1. 

The four-component sum of the photon emission from the test-particle result for SN1006 is 
shown as the heavy dotted curve in Figure ^ where the emission volume has been adjusted to fit the 
radio flux. The broad-band fit is unsatisfactory and does not improve substantially if parameters 
other than ?7inj,p are varied. In the test-particle result, the X-ray emission is dominated by thermal 
bremsstrahlung and the TeV to radio flux ratio is much smaller than in the nonlinear result. 



10 



Of course, many examples exist of full nonlinear calculations of SNR evolution (e.g., Dorfi & Volk 1996; Berezhko 
1999b; Kang & Jones 1995), but despite this, we believe the importance of nonlinear effects has been slow to 



et al. 

reach the general astrophysical community. 
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4. SUMMARY AND CONCLUSIONS 

We have developed a computationally fast and easy-to-use model of nonlinear diffusive shock 
acceleration along with the accompanying photon emission from the resultant electron and ion 
spectra. Using this model, we showed how the emission depends on both model and environmental 
parameters and identified ambient density and magnetic field as being the most important in deter- 
mining the broad-band emission. We also showed that our simple model is in excellent agreement 
with a more complete and complex model of the supernova remnant SN1006. 

A particularly important aspect of nonlinear acceleration is that shock heating is linked to 
particle acceleration and thus to the broad-band photon emission. As shown in Figures g ||a, 
and [ll], shocks which accelerate particles efficiently have lower postshock temperatures and higher 
postshock densities than test-particle predictions. Nonlinear shocks also produce electron spectra 
with superthermal tails (Figures ||, ||, and |TI|). These factors are important for X-ray line emission 
(e.g., Dorfi & Bohringer |1993| ) and open up the possibility of using radio, X-ray continuum, and 
7-ray observations to constrain X-ray line models and vice versa. We believe the model we present 
here is a first step in this process. 

To describe SNRs, we use a model of an expanding, spherical shock wave (Truelove &: Mc- 
Kee 1999| ) to obtain shock parameters as a function of SNR parameters and time. Using these 



parameters, we calculate nonlinear particle soectra and produce "snapshot" continuum photon 
spectra from synchrotron, bremsstrahlung, inverse-Compton, and pion-decay, spanning the range 
from radio to TeV 7-rays. More realistic models would include inhomogeneous ejecta, emission from 
dense clumps, pre-SN winds, a reverse shock, volume integrated emission from regions undergoing 
adiabatic cooling, oblique shock geometry, and effects on the SNR evolution from particle escape 
(or dilution). Some of these additions have already been performed by Berezhko, Ksenofontov, &; 
Petukhov ( |1999a ) and our results are in excellent agreement with theirs for SN1006 (Figure |9| and 



Section |3.2.l[) , some can be added easily to the simple model, but others (such as oblique geome- 
try) are much more difficult. Despite this, we believe this nonlinear model is accurate enough to 
approximate expected emission fluxes and clearly track important trends as we describe in detail 
in Section |3.1.2| . 

Nonlinear shock acceleration unavoidably involves a large number of model and environmental 
parameters and it is essential to know how they relate to each other and which ones have the 
greatest impact on the results. Of the model parameters, the injection efficiency, f7inj, P , has the 
greatest influence on the solutions since it sets the overall efficiency and determines whether the 
acceleration is nonlinear or can be treated as test particle (Figures ^ and ||a). Unfortunately, ??mj,p 
is not well constrained by theory and heliospheric observations remain limited; in fact, there's 
virtually no constraining information on the efficiency of high Mach number shocks like those 
expected in young SNRs. However, the differences between test-particle and nonlinear predictions 
are so large in both the broad-band continuum and the X-ray line emission, that we believe it 
will be possible to set strong constraints on the acceleration efficiency once nonlinear models of 
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X-ray line emission become available. Other model parameters, such as the ratios of shocked 
electron to proton temperature and of electron to proton acceleration efficiency, are less fundamental 
but produce characteristic differences in the emission which can be constrained with sufficiently 
complete models and observations (Figures |I]b and c) . 

The environmental parameters include the ambient interstellar proton density n p Q and density 
profile, the magnetic field Bq, the ejecta mass and density profile, supernova explosion energy, etc. 
Of these, the ambient density and magnetic field have, by far, the largest effect on the broad-band 
emission. As we show in Figures ^, || and [j| varying n p Q and Bq produce a complicated set of 
changes in the emission which cannot be simply characterized. If density is held constant, the 

3 /2 

intensity in the radio band scales roughly as Bq , but at photon energies E > MeV, the intensity 
can decrease with increasing Bq due to nonlinear effects (Figure ||b) and a weakening of the shock 
(i.e., decrease in the Alfven Mach number). Likewise, when Bq is fixed, an increase in density 
generally causes an increase in emission in the MeV range, but the emission in the radio band 
is relatively insensitive to density (Figure ||a). At TeV energies, the maximum photon energy 
increases as n p o is decreased and, in general, the highest photon energies are obtained with the 
largest Bq and the lowest n p Q. The TeV/radio flux ratio, however, is greatest for low ambient 
magnetic fields (Figure [?]) and this is an important parameter if radio SNRs are selected for hard 
7-ray searches. As we show in Figure ||, the distinctive pion-decay bump will be most prominent 
for remnants in high n p Q and high Bq environments; given the flux levels predicted, we anticipate 
positive detections of such spectral features (generally below EGRET sensitivities for n <^ 10 cm -3 ), 
by the GLAST experiment in the not too distant future. While TeV experiments cannot probe 
such spectral features, they can detect the upper range of pion decay emission in remnants in high 
Bq surroundings; this situation may already be realized in the recent announcement (Volk et al. 
2000) of a positive detection of Cas A by HEGRA. 

The nature of the X-ray emission also depends importantly on density and magnetic field. As 
density is increased at a given Bq, the keV X-ray emission goes from being totally dominated by 
synchrotron from relativistic electrons to quasi-thermal bremsstrahlung emission. Note that the 
bremsstrahlung dominates even in the radio band for n p o = 10 cm" - 3 and Bq = 3 //G(Fig.ga). bim- 
ilarly, at a given density (Figure |5|b), the X-ray emission goes from being thermal bremsstrahlung 
to synchrotron as -Bo increases. If synchrotron is dominant, X-ray lines will be weak or absent so 
these differences are readily distinguishable observationally. 

Finally, we suggest that the most important aspect of modeling photon emission from astro- 
physical shocks depends on whether or not the acceleration is efficient and nonlinear or inefficient 
and test-particle in nature. These situations can be quite different and if the acceleration is non- 
linear, all parameters are interconnected and observations in any energy band limit the freedom to 
vary parameters in all other bands. Since there are a large number of parameters, understanding 
the nonlinear interactions requires a model that can efficiently map parameter space. We believe 
the model presented here, while not a replacement for more complete models, can do this expedi- 
ently and accurately and be an aid in interpreting the vast amount of information expected from 
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current and future space and ground-based telescopes. 
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TABLE 1 
Parameters for Figure |TJ and SN1006 



Input parameters 


Low 


High 


oJN 10U6nl 


T>T^TDnn5 
BK.F99 


bJN 1U06tp 


n p0 [cm -3 ] 1 


0.01 


10 


0.1 


0.1 


0.1 


Bo \fiG] 


3 


30 


9 


9 


9 


B 2 [/iG] 2 


3 


30 


30 


30 


30 


*snr [yr] 


1000 


1000 


990 


990 


990 


E SIl [10 51 erg] 


1 


1 


1 


1 


1 


M ej [Mq] 


1 


•5 


1.4 


1.4 


1.4 


T p0 [K] 


10 4 


10 4 


10 4 


10 4 


10 4 




1 


1 


1 




1 


10~ 3 


10~ 3 


2.1 xlO" 4 


5 xlO" 4 


lxl0~ 5 


(e/p)rel 


0.03 


0.03 


2xl0~ 3 


2xl0~ 3 


2xl0~ 3 


?7mfp 


1 


1 


13 




13 


Q: 


1 


1 


0.6 


0.5 


0.6 


n 








7 


7 


7 


Output values 


M so 


400 


85 


235 


170 


235 




120 


77 


72 


67 


72 


V; k [km s _1 ] 


6400 


1350 


3750 


3520 


3750 


-Rsk [pc] 


11 


3.0 


7.0 


6.9 


7.0 


r to t 


9.7 


8.5 


7.2 


6.7 


4.1 


'"sub 


3.0 


2.9 


3.6 


3.6 


4.0 


^max,p [TeV] 


30 


16 


9 


10 


11 


£max,e [TeV] 


30 


12 


9 


10 


10 


^/inj,e 


5xl0" 3 


0.11 


2.9xl0~ 5 




8xl0~ 7 


^ p2 [K] 


6.4 xlO 7 


3.6 xlO 6 


4.7xl0 7 




1.7xl0 8 


Tt P [K] 4 


5.0 xlO 8 


2.2 xlO 7 


1.7xl0 8 




1.7xl0 8 


£rel 


0.80 


0.75 


0.63 


0.40 


0.02 


Flux parameters 


D snr [kpc] 


1 


1 


1.7 


1.7 


1.7 


Vemis [PC 3 ] 


1 


1 


195 6 




5 xlO 3 7 



1 This is the unshocked proton density. In all cases, including our match to the Berezhko, 
Ksenofontov, & Petukhov example, a helium number density of 0.1 n p o is assumed. 

2 In order to compare our model to that of Berezhko, Ksenofontov, & Petukhov ( 1999a| ), 
we include a compressed downstream magnetic field for purposes of calculating the synchrotron 
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emission from SN1006. In all other results we take B2 = Bq. 

3 See the discussion in Section 3.2.1 concerning this parameter and the Berezhko, Ksenofontov, 
& Petukhov results. 

4 This is the temperature of the shocked gas in the test-particle approximation. 

5 The parameters given here are for the model shown as a solid line in Figure Id of Berezhko, 
Ksenofontov, & Petukhov ( |1999a| ). 

6 This value is that given in Eqn. (|19|). 

7 This value is ~ 15 times V em i s given in Eqn. (|l9|) . 
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TABLE 2 
Parameters for Figure || 







(a) 








(b) 






Input 


solid 


dashed 


dot-dashed 


dotted 


solid 


dashed 


dot-dashed 


dotted 


n p0 [cm -3 ] 


0.01 


0.1 


1 


10 


1 


1 


1 


1 


B \pG] 


3 


3 


3 


3 


3 


10 


30 


100 


isnr [yr] 


1000 


1000 


1000 


1000 


1000 


1000 


1000 


1000 


E sn [10 51 erg] 


1 


1 


1 


1 


1 


1 


1 


1 


M ej [M ] 


1 


1 


1 


1 


1 


1 


1 


1 


T p0 [K] 


10 4 


10 4 


10 4 


10 4 


10 4 


10 4 


10 4 


10 4 




1 


1 


1 


1 


1 


1 


1 


1 




10~ 3 


10~ 3 


10~ 3 


10~ 3 


10" 3 


10~ 3 


10~ 3 


10~ 3 


(e/p)rel 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


Output 


M so 


400 


220 


130 


79 


130 


130 


130 


130 


M ao 


110 


200 


370 


720 


370 


110 


37 


11 


V sk [km s -1 ] 


6340 


3440 


2050 


1260 


2050 


2050 


2050 


2050 


#sk [pc] 


10.5 


7.4 


4.8 


3.1 


4.8 


4.8 


4.8 


4.8 


not 


9.6 


12 


14 


17 


14 


9.6 


6.7 


4.8 


'"sub 


3.0 


3.0 


3.0 


3.0 


3.0 


2.9 


3.0 


3.2 


-Emax,p [TeV] 


32 


8 


2.4 


0.8 


2.4 


11 


44 


180 


-Emax,e [TeV] 


32 


8 


2.4 


0.8 


2.4 


11 


20 


12 


^?inj,e 


5xl0~ 3 


0.035 


6 xlO" 3 


7 xlO" 3 


6 x 10~ 3 


9 xlO" 3 


0.011 


0.013 


^p 2 [K] 


6.3 xlO 7 


1.3 xlO 7 


3.0 xlO 6 


7.8 xlO 5 


3.0 xlO 6 


6.6 xlO 6 


1.3 xlO 7 


2.8 xlO 7 


Zt P [K] 


4.9 xlO 8 


1.5xl0 8 


5.1 xlO 7 


1.9 xlO 7 


5.1 xlO 7 


5.1 xlO 7 


5.1 xlO 7 


5.1 xlO 7 




0.81 


0.85 


0.88 


0.90 


0.88 


0.79 


0.63 


0.37 


Flux 


L> snr [kpc] 


1 


1 


1 


1 


1 


1 


1 


1 


Kmis [PC 3 ] 


510 


140 


33 


7.2 


33 


50 


70 


100 
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10" 3 0.01 0.1 1 10 100 1000 10 4 10 5 
Momentum, p/(m c) 

Fig. 1. — Downstream phase space distribution functions, /, versus momentum, p. We have mul- 
tiplied f(p) by [p/ (m p c)] 4 to flatten the spectra, and by [(m p c) 3 /n p o] to make them dimensionless. 
The solid histogram is a Monte Carlo model result, while the dashed, dot-dashed, and dotted curves 
are from the simple model. The three simple model results are obtained with identical input pa- 
rameters except for ?7mj,p which is varied as shown. The ??mj,p = 5x 10 -3 case is chosen to match 
the Monte Carlo result and a good correspondence between the two models is obtained above the 
injection momentum, j>; n j ~ 7xl0~ 3 m p c. The injection momentum, p- m j, varies with input pa- 
rameters: for ??mj,p = 3xl0~ 4 , pi n j ~ 1.0xl0~ 2 m p c, and for r/i n j,p = 5xl0~ 5 , p- m j ~ 0.11 m p c. 
As explained in Berezhko & Ellison (1999), the model spectra are considerably steeper than the 
test-particle prediction for r to t = 46, / oc p~ 3 - 06 (solid line), at the highest momenta. 



-6-3 3 6 

Log 10 Momentum, p/(m c) 




-6 -3 3 6 9 

Log 1Q Energy, [MeV] 



Fig. 2. — Top panel: Phase space distribution functions, f(p), flattened and made dimensionless 
as in Fig. |. The solid curve shows protons, the dashed curve shows electrons with (e/p) ie \ = 0.05, 
and the dotted curve shows electrons with (e/p) re \ = 0.01. No synchrotron or inverse-Compton 
losses are included so the electron and proton spectra cutoff at the same momentum. Bottom 
panel: Differential energy flux distributions for the same particles. These spectral are normalized 
with a far upstream number density, n' , such that u'qUq = 1 cm -2 s~ x . The dot-dashed curves are 
test-particle electron results with (e/p) re \ = 0.01 included for comparison. All input parameters 
are the same as for the dotted curves except ?7i n j,p = 5x 10~ 6 which gives r to t = 4.1. 
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Fig. 3. — Top panel: Phase space distributions for a parameter set with low values of n p o, M e j, and 
Bq (light-weight curves) and high values (heavy-weight curves). The protons are shown as solid 
curves and the electrons are shown as dashed curves. Bottom panel: Photon fluxes at Earth from 
the distributions in the top panel assuming an emission volume of V cm \ s = 1 pc 3 at a distance of 
1 kpc. In each case, the synchrotron emission is shown with a dashed curve, the bremsstrahlung 
with a dot-dashed curve, the inverse-Compton with a dotted curve, and the pion-decay with a solid 
curve. The inverse-Compton emission for the high value example is insignificant and is not shown. 
The full set of parameters for these models is given in Table 1 under columns marked 'High' and 
'Low.' 
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Fig. 4. — Total photon emission for various model parameters. In (a), ?7inj,p = lxlCT 5 (solid), 
lxlO" 4 (dashed), lxlO" 3 (dot-dashed), and 2xl(T 3 (dotted) all with T c2 /T p2 = 1 and (e/p) re i = 0.03. 
In (b): T c2 /T p2 = 0.1 (solid), 0.3 (dashed), and 1 (dot-dashed) all with r? inj)P = 1 x 10~ 3 and 
(e/p)rci = 0.03. In (c): (e/p) rd = 0.01 (solid), 0.03 (dashed), and 0.1 (dot-dashed) all with 
7?inj,p = 1 x 10~ 3 and T C 2/T P 2 = 1. In all cases, n p o = 1 cm" 3 , Bq = 3 [iG, E sn = 1 x 10 51 erg, 
M c j = IMq, and t snr = 10 3 yr. Note that the vertical axis is proportional to the standard v F v 
representation. 
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Fig. 5. — Total photon emission for various input parameters. In (a), the ISM field is fixed at 
Bq = 3 /jlG and the ambient number density is varied such that: n p o = 0.01 cm -3 (solid), 0.1 cm" 3 
(dashed), 1 cm -3 (dot-dashed), and 10 cm -3 (dotted). In (b), Bq is varied: Bq = 3 /j,G (solid), 
10 fiG (dashed), 30 fiG (dot-dashed), and 100 ^G (dotted), with the density pinned to n p o = 1 
cm -3 . The input and output parameters for these examples are given in Table 2, but in all cases, 
r/inj,p = 10~ 3 , T c2 /T p2 = 1, (e/p) rel = 0.03, E sn = 10 51 erg, M ej = 1M , t snr = 1000 yr, a = 1, 
^mfp = 1) an d the flux at Earth is determined assuming D snr = 1 kpc and an emission volume given 
by Eqn. (|l9|). Also depicted are the canonical integral flux sensitivities for Veritas and GLAST 
(e.g., Weekes et al. 1999). 
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Fig. 6. — Total photon emission for various input parameters. In (a), E sn = 0.2 xlO 51 erg (solid), 
1 x 10 51 erg (dashed), and 5xl0 51 erg (dot-dashed). In (b), i snr = 30 yr (solid), 100 yr (dashed), 



10 yr (dot-dashed), and 10 yr (dotted). In all cases, 



10~ 



n 



1 cm" 3 , B = 3/iG, 



T e 2/T p 2 = 1, (e/p) re i = 0.03, M e j = IMq, a = 1, r/ m f p = 1, and the flux at Earth is determined 
assuming D snr = 1 kpc and an emission volume given by Eqn. (W{ 



- 36 - 



T 1 1 1 — I I I I | 1 1 1 1 — I I I I 




1 10 100 

Magnetic Field, B [/xG] 

Fig. 7. — The ratio of gamma-ray to radio fluxes, -Rphotj as denned in the text as a function 
of magnetic field and unshocked density as marked in cm -3 . All plots have E sn = 10 51 erg, 
M ej = 1M , (e/p) rel = 0.03, Tfej.p = 10" 3 , t sni = 1000 yr, a = 1, ry mfp = 1, and T c2 /T p2 = 1. The 
dots indicate parameters used in Figure p| 
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Fig. 8. — Inverse-Compton and pion-decay components for the densities and magnetic fields marked 
with dots in Fig. [7|. 
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Fig. 9. — Photon spectra for SN1006. All spectra are calculated at the shock and show a 'snapshot' 
of the remnant at 990 yr after the explosion. The heavy-weight solid and dotted curves are sums of 
the four photon components. The component spectra for the heavy-weight solid curve are shown 
and labeled. The heavy solid curve is the result for a nonlinear shock, while the heavy dotted 
curve is a test-particle result. The radio data are from Reynolds & Ellison ( |1992 ) and the X-ray 
data (solid points) are adapted from Reynolds ( |1996 ), the EGRET upper limit (cross) is from 
Mastichiadis & de Jager (1996), and the CANGAROO TeV points (squares) are from Tanimori 
et al.(1998). Note that the integral CANGAROO points are plotted, with no adjustment, on our 
differential representation. 
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Fig. 10. — Proton phase space distributions from the simple model calculated here (dashed line) 
and the kinetic model of Berezhko, Ksenofontov, & Petukhov (1999a) (solid line) for SN1006. The 
Berezhko et al. result is the solid curve in their Figure Id, and no adjustment in normalization is 
made to improve the match. 
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Fig. 11. — Electron and proton phase space distributions for the nonlinear and test-particle models 
shown in Fig. |^. The light-weight curves show the nonlinear f(p), while the heavy-weight curves 
show the test-particle f(p). All input parameters, other than the proton injection efficiency, are 
the same in the two cases. Compared to the nonlinear result, the test-particle example has a higher 
shocked temperature (as indicated by the position of the thermal peak) and a lower shocked density. 
The difference in density is not apparent from the figure because of the normalization. 



